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ABSTRACT 

The afterglow of GRB 050401 presents several novel and interesting features : 

(i) An initially faster decay in optical band than in X-rays. 

(ii) A break in the X-ray light curve after ~ 0.06 day with an unusual slope after 
the break. 

(iii) The X-ray afterglow does not show any spectral evolution across the break 
while the R band light curve does not show any break. 

We have modeled the observed multi-band evolution of the afterglow of GRB 
050401 as originating in a two component jet, interpreting the break in X-ray light 
curve as due to lateral expansion of a narrow collimated outflow which dominates 
the X-ray emission. The optical emission is attributed to a wider jet component. Our 
model reproduces all the observed features of multi-band afterglow of GRB 050401. 

We present optical observations of GRB 050401 using the 104-cm Sampurnanand 
Telescope at ARIES, Nainital. Results of the analysis of multi-band data are presented 
and compared with GRB 030329, the first reported case of double jet. 
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1 INTRODUCTION 

The optical and X-ray light curves of Gamma Ray Burst 
(GRB) afterglows, in the simplest cases, show a power law 
decay with an index a ~ 1.0. Deceleration of the relativis- 
tic shock wave generated by the explosion which results in 
GRB can explain the power law decay of the GRB after- 
glows. The most common deviation from the power law de- 
cay behaviour of the afterglow light curves is an achromatic 
break seen in the light curve. This break has been seen in a 
significant number of GRB afterglows and has been success- 
fully explained as being due to the sideways expansion of 
the collimated ejecta from the explosion. In the post Swift 
era, many more deviations from this simple behaviour of the 
afterglow light curve have been detected. Swift with its ca- 
pabilities of quick slewing towards the source has been able 
to observe GRB afterglows as early as a few tens of seconds 
after the burst. In this early part of the evolution the GRB 
afterglows commonly exhibit a steep decay with a ~ 3 to 
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5 with the usual definition F v (t) oc t 



where F„(t) is 



the observed afterglow flux at frequency 'u' and time t. The 
phase of steep decay lasts for about a few hundred seconds 
after which a slower decay, with a ~ 0.5, of the afterglow 
starts. About a few thousands of seconds after the burst the 
afterglow starts decaying steeply again with a > 1.0. 

Many GRB afterglows observed by Swift show puz- 
zling features in the light curves like (1.) early steep de- 
cay [a ~ 3 to 5] and (2.) Chromatic breaks (breaks seen 
in some wavebands but not others) with Aa ~ 1.0 which 
are difficult to explain using the standard firebal l model 
(|Rees and Meszaroslll992l; iMeszaros and Reeslll993l). It has 
been shown bv lO'Brien et all l|2006l ); IWillingale et al.l(|2006h 
that the puzzling features of the X-ray afterglow light curves 
can be fitted using one or two components with exactly the 
same empirical functional form, viz. an exponential fall fol- 
lowed by a powerlaw decay of flux with time, although it 
has not yet resulted into any physical understanding of the 
behaviour of the X-ray afterglow. While there is no clear 
understanding of the early steep decays of GRB afterglows, 
a few plausible explanations have been put forward : see 
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e.g. IZhang et al] (|2006l ); iPe'er et all (|2006l ). The flat decay 
of X-ray afterglow light curves which follows the steep de- 
cay have been, in some cases, explained as being due to en- 
ergy injection from the central en gine, probably a magnetar 
llZhang and MeszaroslboOll . 12002! ). From the study of chro- 
matic breaks seen in six well sampled afterglow light curves 
IPanaitescu et all l|2006l ) concludes that if both, the optical 
and the X-ray afterglows, were to arise from the same out- 
flow then the chromaticity of light curve breaks can rule out 
energy injection or the structure of the jet as the possible 
reasons of it. 

One such GRB afterglow with puzzling features in 
optical and X-ray light curves is GRB 050401. GRB 
050401 tr iggered Swift - B AT at 14:20:15 UT on 2005 
April 01 i|Barbier et al.l 120051). The X- r ay af terglow was 
detected by Swift-XRT i|Angelini et al.1 120051 ) about 130 
seconds after the trigger and the optical afterglow can- 
didate was confirme d by g round based observations by 
iPrice and McNaughtl l|2005h . The burst duration Tgo is 
estimated to be ~ 33 seconds (|Sakamoto et al.ll2005l ). Using 
the measured spect ral redshift of the afte rglow (z = 2.9) 
dFvnbq et all 120051) and the fluence (ISakamoto et all 
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the isotropic equivalent energy released during the explo- 
sion turns out to be 1.4 x 10 54 for a flat universe with 
fl m = 0.3, S1a = 0.7 and H = 70 km s" 1 Mpc" 1 . 

Multiband afterglow of GRB 050401 also presents some 
puzzling features which can be summarized as follows : 

(i) A break in the X-ray light cu rve after ~ 0.06 day with 
an unusual slope a fter the break (|De Pasauale et alj|2006l : 
IWatson et a"ill2005l ). 

(ii) The X-ray afterglow does not show any spectral evo- 
lution across the br eak while the R band light curve does 
not s how any break (|De Pasauale et al]|2006l : IWatson et al.l 
120051 ). 

(iii) A large extinction inferred from X-ray afterglow 
which is not consist ent with the observed optical afterglow 
jWatson et al.ll2005l) . 

The optical observations are presented in § [2] We have 
done some preliminary analysis of the light curves which is 
discussed in § [3] We have tried to explain the multi-band 
behaviour of the GRB afterglow using a double jet model 
which is described in § U along with the previous attempts 
by others using a different model. In the Discussion section 
(§ [5| molecular clouds as a plausible explanation for the 
large extinction is presented (§ 15. ip . The only other GRB 
afterglow which has been ex plained using a simi l ar double jet 
mode l is the GRB 030329 (|Berger et al.ll2003l : iResmi et all 
l2005l ). We compare the physical features of GRB 030329 
and GRB 050401 in § 15.21 Our conclusions are summarized 
in § El 



2 OPTICAL OBSERVATIONS AND DATA 
REDUCTION 

Optical observations of the afterglow of GRB 050401 were 
carried out in the broad and Johnson V and Cousins RI fil- 
ters using the 104-cm Sampurnanand Telescope of ARIES, 
Nainital on 01 April 2005. The gain and read out noise of the 



CCD camera are 10 e~/ADU and 5.3 e~ respectively. The 
data have been binned in 2x2 pixel 2 to improve the signal- 
to-noise ratio. The bias subtracted, flat fielded and cosmic 
ray removed images were process ed and analyse d using MI- 
DAS 0, I RAFBan d DAOPHOT (|Stetsonll 19871 ) softwares. 

The lLandoltl (|l992l ) standard region SA 107 and the 
OA field in BVRI filters was observed on 16 May 2005 for 
photometric calibration during good photometric sky con- 
ditions. The values of atmospheric extinction on the night 
of 16/17 May 2005 determined from the observations of SA 
107 bright stars are 0.26, 0.18, 0.13 and 0.10 magnitude in 
B, V, R and I filters respectively. The 7 standard stars in the 
SA 107 region cover a range of 0.339 < (V - R) < 0.923 in 
color and 12.116 < V < 14.884 in brightness. 

Using these transformation coefficients we determine 
BVRI magnitudes of 18 secondary stars in GRB 050401 field 
and their average values are listed in Table [1] The (X, Y) 
CCD pixel coordinates were converted to 02000, £2000 v alues 
using the astrometric positions given bv lHendenll|2005l ). The 
18 secondary stars in the field of GRB 050401 were observed 
2 to 4 times in B, V, R and I filters. These stars have internal 
photometric accuracy better than 0.01 mag. The zero-point 
di fferences on co mparison between our photometry and that 
of lHendenl (|2005l ) are 0.15 ± 0.08, 0.09 ± 0.04, 0.10 ± 0.05 
and 0.54 ± 0.29 magnitude in B,V,R and I filters respec- 
tively. These differences are based on the comparison of the 
6 secondary stars in the GRB 050401 field. 

The afterglow magnitudes were differentially calibrated 
with respect to the secondary stars listed in Table [1] The 
magnitudes derived in this way are given in Table [5] 



3 LIGHT CURVES OF GRB 050401 
AFTERGLOW 

Along with our own observations we have used observa- 
tions reported elsewhere to study the light curves of GRB 
050401. Th e X-ray light curve of GRB 050401 was ob- 
tai ned from Watson et aTl l|2005l ). The optical observations 
by I Watson et all (|2005l ) have been calibrated by observ- 
ing a Landolt field. We do not have a detailed informa- 
tion about this calibration. Hence, to take into account any 
uncertainties associated with it, we have added an error 
of 0.2 magnitudes i n the optical observations reported by 
I Watson et aL ( 20051). Another set of optical observations is 
taken from Rvk off et al. I (|2005l ) the calibration of which is 
roughly equivalent to the R c band system. We add a small 
error of magnitude 0.1 to all these observations by ROTSE- 
III to take into account the calibration uncertainties. 

VLA report ed a 4<t detection of a source at the position 
of GRB 050401 l|Soderberdl2005l ) with intensity of 122 ^Jy 
at 8.46 GHz about 5.7 days after the burs t. Other attempts, 
inclu ding by GMRT in India at 610 MHz l|Chandra and Ravi 
120051) and by ATCA in Australia at 8.5 GHz and 4.8 GHz 
i Saripalli etai1l2005h . to observe the radio afterglow of GRB 
050401 could produce only upper limits. 

To construct the optical light curve we have corrected 



1 MIDAS is distributed by the European Southern Observatories. 
Visit : www.eso.org/esomidas/ 

2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, USA. Visit : |http:/ /iraf.noao.edu/| 
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_steeper after the break. 



ID 



C«2000 
(h m s) 



<52000 
(deg m s) 



V 
(mag) 



B — V 
(mag) 



V -R 

(mag) 



V — I 
(mag) 



1 


16 


31 


20.01 


02 


06 


52.9 


17.28 


0.64 


0.35 


0.84 


2 


16 


31 


23.84 


02 


07 


44.3 


16.87 


0.77 


0.52 


0.98 


3 


16 


31 


29.22 


02 


08 


13.8 


17.66 


1.14 


0.75 


1.40 


1 


16 


31 


37.63 


02 


08 


07.3 


16.78 


0.69 


0.46 


0.87 


5 


16 


31 


40.12 


02 


10 


30.1 


16.35 


0.58 


0.38 


0.73 


6 


16 


31 


36.96 


02 


11 


36.5 


18.23 


0.97 


0.62 


1.12 


7 


16 


31 


32.61 


02 


12 


38.7 


17.68 


0.41 


0.34 


0.67 


8 


16 


31 


24.79 


02 


13 


35.4 


19.56 


1.37 


1.08 


2.29 


9 


16 


31 


18.94 


02 


13 


12.1 


19.21 


1.26 


0.83 


1.60 


10 


16 


31 


18.56 


02 


12 


40.8 


15.30 


0.85 


0.51 


0.93 


11 


16 


31 


22.46 


02 


11 


13.7 


15.61 


0.69 


0.46 


0.87 


12 


16 


31 


21.38 


02 


10 


43.0 


15.51 


0.88 


0.53 


0.99 


13 


16 


31 


19.42 


02 


09 


56.3 


14.60 


0.55 


0.36 


0.69 


11 


16 


31 


15.08 


02 


09 


19.1 


14.34 


0.61 


0.38 


0.74 


15 


16 


31 


23.42 


02 


09 


13.7 


16.39 


0.66 


0.44 


0.83 


16 


16 


31 


17.26 


02 


07 


58.9 


16.17 


0.63 


0.41 


0.81 


17 


16 


31 


15.93 


02 


07 


36.6 


18.91 


1.16 


0.89 


1.75 


18 


16 


31 


14.79 


02 


07 


14.6 


17.54 


0.74 


0.47 


0.93 



Table 1. The identification number (ID), (a, S) for epoch 2000, 
standard V, (B — V), (V — R) and (V — /) photometric magnitudes 
of the stars in the GRB 050401 region are given. 



Date (UT) 
2005 April 


AT 

(days) 


Magnitude 
(mag) 


Passband 


01.8824 


0.2850 


22.33 ± 0.347 


V 


01.8324 


0.2850 


21.43 ± 0.231 


R 


01.8698 


0.2724 


20.51 ± 0.207 


I 



The X-ray and optical (R band) afterglow of GRB 
■050401 is very well sampled over a wide period of observa- 
tion. The available R band observations cover a duration 
from 36 s to 13 days after the burst while the X-ray 
observations range from ~ 130 s to 12 days after the burst. 
The X-ray afterglow light curve shows a prominent break 
near 0.06 day while the optical afterglow does not show any 
such break in the light curve. We analyse this behaviour in 
detail below. 



1. The X-ray light curve shows a clear break near 0.06 day. 
The change of slope across the break is significant. Fitting 
Equation [T] to the data yield the decay slopes 

Qxi = 0.58 ± 0.02 for At < 0.06 day; 
ax2 = 1.37 ± 0.03 for At > 0.06 day; 

The change of slope across the break, A«x ~ 0.8, is there- 
fore quite substantial. 

2. The optical afterglow shows a monotonic decay with 
_decay index «r = 0.82 ± 0.02 over the entire period of ob- 
servation (up to 13 days). There is no evidence of a break 
simultaneous with that in the X-ray light curve. 

3. According to the standard fireball model of GRB after- 
glows the X-ray light curve is expected to decay at least as 
fast as the optical light curve which is indeed true for the 
majority of GRB afterglows observed so far. In the case of 
GRB 050401, we find that the X-ray afterglow shows a decay 
slower than optical light curve till ~ 0.06 day after which it 
decays at a much faster rate as described above. 

Thus, the relatively slow initial decay of optical and X- 
ray light curves, presence of a break in X-ray light curve and 
absence of such a break in optical, and initially slower decay 
of the light curve in X-rays than in optical bands makes the 
afterglow of GRB 050401 an unusual and interesting one. 



Table 2. The optical observations of the afterglow of GRB 050401 
using the 104-cm Sampurnanand Telescope at ARIES, Nainital. 
AT in column 2 refers to the time after the burst in days. The 
effective exposure time after combining all the images turns out 
to be 900 s for individual passbands reported here. 

the observed ma gnitudes for th e standard Galactic extinc- 
tion law given bv lMathisl(|l990h . The Galactic extinction in 
the direction of GRB 050401 is estimated to be E(B-V) = 
0.065 mag from the smoothed reddening map provided by 
ISchlegel et al.l (1 19981). The effe c tive w avelength and normal- 
isation given by Bessell et al.l (1 19981 ) were used to convert 
the magnitudes to fluxes in pjy. 

Most of the GRB afterglow light curves are well char- 
acterised by a broken power law of the form 

F = F {(t/t b )^ 3 + {t/t b T 23 y 1/s (i) 

where «i and 0:2 are the afterglow flux decay indices before 
and after the break time (tb), respectively. Fo is the flux 
normalisation and 's' is a smoothening parameter which 
controls the sharpness of the break. Most known GRB 
afterglows have ai ~ 1 and ai > a.\ i.e. the decay becomes 



4 MODELING OF GRB 050401 AFTERGLOW 

The change in slope across the break in the X-ray light curve 
Aqx ~ 0.8 is too large to be explained by the passage of 
a spectral break. In the standard fireball model of GRB 
afterglows the passage of the cooling break v c through the 
observing band leads to a steepening of light curve by an 
amount A«x = 0.25, much smaller than that is observed 
for GRB 050401 afterglow, along with the change of spectral 
slope by A/3 X = 0.5. The X-ray spectrum of GRB 050401 
does not exhibit any change in the spectral slope across the 
break. We thus rule out the possibility of v c passing through 
the X-ray band at the t ime o f break. 

|Pe Pasquale et aLI (|2006l ) explains the initial flatter de- 
cay and the break in the X-ray light cur ve based on a 
model bv IZhang and Meszarosl l)200ll . [2002). According to 
this model, the central engine of GRB remains active for sev- 
eral thousand seconds after the burst, continuously inject- 
ing energy into the fireball. If the central engine is injecting 
energy above a certain critical rate then it can slow down 
deceleration of the shock wave which results in a shallow 
decay of the light curve. The break in the light curve occurs 
when the central engine stops injecting sufficient amount of 
energy into the fireball. After this epoch the afterglow can be 
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described using standard fireball model giving a = (3/2)/3. 
Being a dynamical effect, the end of energy injection episode 
would result in an achromatic break in the afterglow light 
curves. Although this model seem to explain the X-ray light 
curve reasonably well, absence of a similar break in optical 
afterglow light curve is sufficient to rule this model out for 
GR B 050401. 

IWatson et al" I i|2005l ) point out another puzzle : the soft 
X-ray absorption implies an equivalent optical extinction of 
magnitude A v = 9.11* g magnitudes in the host galaxy, as- 
suming solar abundance. 

However, if the optical and the X-ray emission are part 
of the same synchrotron spectrum, then A v is constrained 
to be ~ 1.45 for no spectral break between optical and 
X-rays and A v < 0.67 if a cooling break exists in be- 
tween (an SMC extinction law is assumed). These values are 
highl y discordant with t hat predicted from X-ray absorp- 
tion. [WatsoiLetjil] (|2005l ) suggests that this may indicate a 
non-universal dust to metals ratio which they estimate to 
be more than a fa c tor of 10 less than that in the SMC. 

IWatson et al.l l|2005h remark that the only alternative to 
this highly anomalous dust to metal ratio is separate emis- 
sion regions for the optical and X-rays. We explore this possi- 
bility assuming that two distinct jet components give rise to 
the observed emission in these two (X-ray and optical) wave- 
length bands. The jet contributing to the X-ray emission is 
narrow, exhibiting an early break while that contributing to 
the optical emission is wider. The optical contribution from 
the narrow jet is strongly diminished due to the presence 
of high extinction A v ~ 9 along the line of sight, while the 
wider jet suffers from a smaller degree of average extinction. 

4.1 Spectral Parameters of the Afterglow of GRB 
050401 

The radiation spectrum of a GRB afterglows exhibits a 
power law spectrum characterised by three break frequen- 
cies - the self absorption frequency v^, the peak frequency Vm 
corresponding to the lower cutoff (7m) in the electron energy 
distribution (n(7<.) oc 7,T P , 7 e > 7 m where y e is the Lorentz 
factor of the radiating electrons) and the synchrotron cooling 
frequency v c . Th e flux F m at z^ m provides the normalisation 
of the spectrum (|Sari et al.lll998l ). 

The photon index (F) of the X-ray afterglow is related 
to its spectral index (/3), F — 1 = /3, which in turn is related to 
the electron energy distribution index p in any given spectral 
regime (f3 = p/2 if v c < ux and (3 = (p — l)/2 if vx < 
Vc). The corresponding temporal decay index ax would be 
(3p — 2)/4 and 3(p — l)/4 respectively before the jet break 
and p in both spectral regimes after the jet break, according 
to the standard fireball model for an afterglow expanding 
in a homogeneous interstellar medium. In the present case, 
the observed values of ax are consistent with p — 1.42 and 
a jet break around 0.06 days after the burst. However, we 
note that the ob served value of the spectral photon index 
T ~ 1.85 ± 0.03 I Watson et al.l 120051 ) implies a steeper p ~ 
1.7 ± 0.06. It should also be noted tha t from analysis of 
the sa me data set of X-ray observations, |Pe Pasquale et ail 
|2006l ) infer (3 = 0.75 ± 0.15 for PC mode data after the 
break at 0.06 days. This /3 is consistent with p = 1.42 that 
we inferred above. 

The optical (R-band) afterglow, on the other hand, 





Narrow Jet 


Wider Jet 


fin (Hz) 

MHz) 

Fpeak(MJy) 

tjet (day) 

p 

E(B - V) H ost 


a-oiS:2i x io 13 

4.1 ± 0.9 x 10 14 
2140+ 210 
0.06 ± 0.03 
1.42 ±0.02 
4.1 


x 10 13 
5.25+ 30 ° x 10 15 

J-' OU -950 

9 -1+0.2 
z - -0.11 

0.23lg;» 




1.2 (85) 



Table 3. Best fit Spectral parameters for the afterglow of GRB 
050401 assuming two component jet model described in §|4] Light 
curves generated using these parameters and their subsequent 
evolution according to the standard fireball model are plotted in 
Figure [l] All the parameters are calculated at 0.01 day after the 
burst. GRB 050401 was at redshift z = 2.9 

exhibits a temporal slope art = 0.82 which, in the com- 
monly encountered spectral regime of i/ m < !/r < v c , im- 
plies p — 2.1. This is different from that inferred for the 
X-ray afterglow, and indeed regardless of the choice of spec- 
tral regimes it is not possible to produce both ax and an 
from the same underlying power-law energy distribution of 
injected electrons. One possibility, therefore, is that the op- 
tical and the X-ray afterglows originate in physically distinct 
outflows. We consider two physically distinct components of 
the outflow, such as the co-axial jets, one having a domi- 
nant contribution in the optical and the other in the X-rays, 
giving rise to the observed afterglow of GRB 050401. 

We then fit the full, multi-band light curves of 
GRB050401 with those predicted by a double jet model us- 
ing linear least square method. Results of this fit are dis- 
played in Figure [T] and the best fit values of various spectral 
parameters are listed in Table [3] We note that the contribu- 
tion of the narrow jet to the optical afterglow is strongly sup- 
pressed due to large extinction. The X-ray afterglow, on the 
other hand, is modified as a sum of the emission from both 
the jets, with the narrow jet being the dominant contributor. 
For the narrow jet we find a best fit value of p — 1.42. For 
the wider jet, which dominates the optical afterglow of GRB 
050401, we estimate p = 2.1. The extinction that the radia- 
tion from the narrow jet encounters is fixed at A v = 9.1 a s 
derived from the soft X-ray absorption ( Wa tson et al.|[2005l ) , 
while that for the wide jet is treated as a fit parameter. 



4.2 Physical Parameters for GRB 050401 

Four spectral parameters (i/ a , v m ,v c and F pca k) are related 
to four physical parameters viz. n (number density of the cir- 
cumburst medium) , E (total energy content of the fireball) , 
energy fraction in relativistic electrons e c and that in mag- 
netic field es. The typical value of self absorption frequency 
i/ a lies in radio-mm waves and hence is best estimated only if 
the afterglow is well observed in these bands. Unfortunately, 
the afterglo w of GRB 05040 1 was detected only once at the 
radio band !|Soderberdl2"005l ) which is not sufficient to deter- 
mine i/ a accurately. We therefore converted the three remain- 
ing spectral parameters into the four physical parameters us- 
ing e e = as an add itional constraint . The choice of this 
relation is motivated bv lMedvedevl (2006). When p < 2.0, as 
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Figure 1. The observed optical & X-ray light curves (left panel) and radio light curves (right panel)of GRB 050401 afterglow compared 
with the double jet model fit (solid lines). The steepening of X-ray afterglow light curves at 0.06 day after the burst is explained as a jet 
break due to the lateral expansion of a narrow jet which has a dominant contribution in X-rays. The surrounding wider jet contributes 
dominantly in optical. Since, no break in the optical light curves is observed till 13 days after the burst, the wider jet is expected to be 
> 29°. Our best fit model gives the value of electron energy distribution index within the narrow jet to be p = 1.42 and that within the 
wider jet to be p = 2.1. The peak in the optical light curves corresponds to the passage of v m through the observing band. The radio 
upper limits are indicated by filled triangles in the right panels. The sole radio detection at 8.5 GHz is indicated by a filled circle. The 
solid lines in the right panels are the light curves expected from the best fit spectral parameters. The corresponding frequencies are listed 
in a rectangle at the top right corner of each box. 



it is in the present case of narrow jet, a high energy cut-off for 
the electron energy distribution is requ ired and the expres- 
sions f or spectral parameters, as given in lWiiers and Galamal 
(1999), have to be mo dified accordin g ly. Th e modifications 
have been provided by iBhattacharval (|200ll ) which we have 
used for estimating the physical parameters in the present 
case. We estimate the density of the circumburst medium 
to be n ~ 10 and e c — y / e~B = 0.03 for both the jets. The 
physical parameters estimated for both the jets are listed in 
Table [4] Using the E lso and n, and the jet break time in 
X-rays, tj — 0.06 days, we find the opening angle of the nar- 
row jet to be quite small, 1.15°. Since there is no jet break 
seen in the optical light curve till ~ 13 days, a lower limit on 
the opening angle of the wider jet is derived to be 29°. The 
collimation corrected kinetic energies are E^ ide > 6.5 x 10 50 
ergs and En arrow — 1.1 x 10 50 ergs. 



Narrow Jet Wider Jet 



n 


14 7+10-5 


zu -19.3 


«e 


(2.3 ±0.6) x 10~ 2 


(4+^) X 10~ 2 




5tt x 10- 4 


lto.9 x !0" 3 




53.23 ± 16.2 


1 Q4+1-36 
1 -°^-0.82 




1.15° ±0.15° 


> 29° 


rpcorr 
^52 


(1.1 ± 0.2) x 10 50 


> 6.5 x 10 50 



Table 4. The physical parameters for the afterglow of GRB 
050401 assuming a two component jet model described in §|4] The 
quantity E^" is the isotropic equivalent energy in units of 10 52 
ergs. The corresponding collimation corrected energy is E^2 rr in 
units of 10 52 ergs. 
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Figure 2. A side view of the double jet (not to scale). The observer is on the axis of the jets and at a distance of 24 Gpc (which can 
be considered as at infinity for geometric purposes in this figure). The arrows indicate the afterglow light rays emanating from the jets. 
The intervening molecular cloud, of size larger than 5 to 55 parsecs, responsible for the observed large extinction is sitting at a distance 
of about 100-1000 parsecs from the GRB. The estimated diameters of the jets around 0.05 days turn out to be about 2 X 10 — 3 pc 
and > 2 X 10 — 2 pc respectively for the narrow and the wide jet. The large cloud covers a significant portion of the central narrow jet 
and partially covers the wide jet when seen from the observer's point of view. As a result, the optical radiation from the narrow jet is 
completely extincted. Most of the optical radiation from the wide jet does not suffer from this extinction. 



5 DISCUSSION 

5.1 A plausible explanation for the large 

extinction inferred from X-ray absorption 

It is now well established from the observations of GRB hosts 
that long G RBs preferentially occur in m assive star forming 
regions e.g. IWooslev and Blooml (|2006T ). The massive star 
forming regions host large molecular clouds. Typical col- 
umn densities of cold molecular clouds are > 10 22 cm -2 , 
densities 100 — 10 4 cm -3 and sizes ~ 20 pc. Giant molecular 
clouds are even denser (10 4 -10 7 cm ) and larger (~ 100 pc) 
l|Shorell2002l ). It is possible that one such cloud in the host 
galaxy of GRB 050401 happens to fall along our line of sight 
which can explain the large extinction inferred from the X- 
ray spectrum. We consider the possibility of radiation from 
the double jet of GRB 050401 being obscured by a molec- 
ular cloud so aligned that it covers the narrow jet of GRB 
050401 completely while the wide jet is partially covered. 
By changing the fractional coverage of wide jet by the cloud 
we measured change in the value of reduced x 2 of the fit. 
In effect, this amounts to adjusting the intrinsic luminos- 
ity of the wide jet upwards with increasing covering factor 
to match the observed optical flux. This results in the rel- 
ative contribution of the wide jet to the X-ray afterglow to 
increase, affecting the fit quality. Keeping all other parame- 
ters fixed at their best fit values obtained for zero coverage, 
we find that a covering fraction of 60% can be accomodated 
within a range of A\ 2 /dof = 1. Beyond this the reduced \ 2 
rises sharply and reaches Ax 2 /dof > 15 for a covering fac- 
tor of ~ 90%. For the observed column density of 1.7 x 10 22 
cm Pasauale et al . 2006), and assuming typical densi- 

ties (100— 1000 cm -3 ) of the molecular clouds, the size of the 
molecular cloud could be estimated to be around 5 — 55 par- 
secs. It is therefore probable that one such molecular cloud 
partially obscures our view of GRB 050401. This situation is 
illustrated in Figure [2] At this point, we would like to point 



out two possible caveats in the double jet model proposed 
here : 

The separation of the optical and X-ray emitting re- 
gions, as proposed in the present model, is motivated by 
the large discrepancy of about 8 magnitudes between the 
amount of optical extinction inferred from soft X-ray ab- 
sorption and that from observed optical-IR spectrum of the 
GRB050 401 afterglow. It shou l d, how ever, be kept in mind 
that the IPredehl and Schmit j (|l995l ) relation used to pre- 
dict Ay from X-ray absorbing column Nh is an empirical 
one, and cannot be considered fully reliable in all circum- 
stances. For example, a metallicity higher than solar by a 
factor of 10, or a dust-to-gas ratio lower by a similar factor, 
can reconcile the X-ray absorption with observed optical ex- 
tinction. Such explanations in this c ase cannot be ruled o ut, 
and have been already discussed bv lWatson et all (|2005l ). 

The second caveat is that the model presented here re- 
quires a rather special geometrical alignment - the two jets 
of the GRB should shine through the outer edge of a molec- 
ular cloud, much larger in size than the transverse extent 
of the jet working surface, in such a manner as to provide 
large extinction to the inner jet but much less to at least 
half the outer component. This requires that the outer edge 
of the cloud be dense, and have a strong density gradient to 
differentially affect the two jet components. An elongated, 
cigar-shaped cloud with its axis nearly parallel to the line 
of sight, would also help such a scenario. We also note that 
the size of the cloud required, as estimated by us using an 
average density, is prone to large uncertainties if its shape 
is unusual or if large density gradients are present. 

5.2 GRB 050401 and GRB 030329 : A comparison 

The only other GRB whose afterglow has been explained as 
being due to double jet is the GRB 030329 (|Berger et al.l 
120031 : iResmi et al.l [200 5) . optical and X-ray light curves of 
GRB 030329 afterglow showed a near simultaneous break at 
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0.55 day whereas the ra dio light curves had a break at about 
10 days after the burst. iBerger et al.l (|2003T ) have explained 
the two breaks as being due to lateral expansion of the two 
co-axial jets of different opening angles (~ 5° and ~ 17°). 

In the case of GRB 050401, afterglow light curves do not 
show the presence of two different breaks. Instead, absence 
of a break at optical frequencies till late times (~ 13 days 
after the burst) leads us to infer the presence of a wider 
jet with opening angle larger than 29° while a steep break 
(Aq ~ 0.8) at 0.06 day after the burst in X-ray light curve 
can be explained as a jet break due to lateral expansion of 
a narrow jet of opening angle 1.15°. 

The wider jet of GRB 030329 was estimated to 
be marginally more energetic tha n the narrower jet 
l|Berger et al.ll2003l ; lResmi et al.ll2005l ). Similarly, in the case 
of GRB 050401, we find, that the wider jet is marginally 
more energetic than the narrower jet. 



5.3 GRB 050401 and the Ghirlanda Relation : 

It has been found that the collimation corrected energies 
(i? 7 ) of the GRBs are correlated with the peak energy of 
the GRB spectrum as measured in the frame of reference 
of the source (E^ k ). This correlation is a lso called as the 
Ghirlanda relation (Gh irlanda et all I2004T ) . Unfortunately, 
the E°lak f° r GRB 050401 is not av ailable as it falls o utside 
However, ISato et ail l|2007h have 

obs 



the energy range of BAT 

used the Konus- Wind spec tral data to find E°^ k . From their 
analysis ISato et ail (|2007T l finds that in order to satisfy the 
Ghirlanda relation the afterglow light curve of GRB 050401 
should exhibit a jet break ~ 10 4 s after the burst. This 
lower limit of the allowed range for jet break time is close to 
the break seen at 0.06 day in the X-ray light curve of GRB 
050401, which we interpret as a jet break corresponding to 
the narrow jet i n our model. 

ISato et al l (|2007l ) quantifies the Ghirlanda relation as 



E 



peak 



a e: 



7,52 



where E. 



7,52 



is the collimation corrected 



energy released in 7 rays during the burst, in u nits of 10 52 
ergs. U sing a sample of a large number of GRBs lSato et all 
( 2007) constrains the value of the proportionality constant A 
: 1950 < A < 4380. Using the estimated value of E\ ao ~ 10 54 
ergs and the 1.15° as the opening angle of the narrow jet in 
our double jet model, the E 1 turns out to be 2 x 10 50 ergs. 



,.„ - U7±HkeV for GRB 050401 as reported by 
ISato et al l pOOTt ) along with _E 7 = 2 x 10 50 ergs we estimate 



Using E. 



src 



A = 70761 llll . This value of A is within 2 a of A = 4380, the 
higher limit on A obtained considering the sample of GRBs 
satisfying the Ghirlanda relation. Having discussed this, we 
would also like to point out that the Ghirlanda relation 
has sometimes been critisized as being due to selection ef- 
fects rather than being an intrinsic correlation l|Butler et al.l 
[200Sh . 



6 SUMMARY 

We have reported VRI band observations of GRB 050401 
afterglow on 1st Apr. 2005. Also, we have modeled the af- 
terglow of GRB 050401 as due to two physically distinct 
collimated outflows, using our own VRI band photometry 
along with the observations available in the literature, and 



compared with GRB 030329. Our main conclusions about 
GRB 050401 are as follows : 

1. We showed that the light curves of GRB 050401 after- 
glow can not be explained under the assumption of contin- 
uous energy injection. The flatter decay, which appealed for 
the continuous energy injection model, can instead be ex- 
plained by low values of electron energy distribution index 
V- 

2. The afterglow of GRB 050401 can be well fit by the 
double jet model with the interpretation that the break in 
the X-ray light curve at ~ 0.06 day after the burst is due to 
a narrow collimated jet expanding sideways. The obscured 
optical emission is attributed to a wider which did not un- 
dergo significant sideways expansion until at least ~ 13 days 
after the burst. 

3. Kinematically, we find that the wider jet is slightly more 
energetic, than the narrow jet. This result is similar to what 
was found in the double jet of GRB 030329. 

4. Our interpretation of the break in the X-ray light curve 
at 0.06 days after the burst as a jet break is consistent with 
the Ghirlanda relation. 
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